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Inventions  (DD882) 


■Modeling  Complex 
Nonlinear  Optical 
Systems 


Alejandro  Aceves 

Department  of  Mathematics  and  Statistics 
The  University  of  New  Mexico 


statement  of  the  problem: 

Study  the  dynamics  of  pulses  in 

photonic  structures. 

What  is  a  photonic  structure  ? 

■  It  is  an  “engineered”  optical  medium  with 
periodic  properties. 

■  Photonic  structures  are  built  to  manipulate 
light  (slow  light,  light  localization,  trap  light) 

■  I’ll  show  some  examples,  but  this  talk  will 
concentrate  on  2-dim  periodic  waveguides 


Motivation  of  study 

We  take  advantage  of  periodic  structures  and  nonlinear  effects  to  propose  new 
stable  and  robust  systems  relevant  to  optical  systems. 


-  Periodic  structure  -  material  with  a  periodically  varying  index  of  refraction  (“grating”) 

-  Nonlinearity  -  dependence  of  the  refractive  index  on  the  intensity  of  the  electric  field 

(Kerr  nonlinearity) 


The  effects  we  seek: 

•  existence  of  stable  localized  solutions  -  solitary  waves,  solitons 

•  short  formation  lengths  of  these  stable  pulses 

•  possibility  to  control  the  pulses  -  speed,  direction  (2D,  3D) 


Prospective  applications: 

•  rerouting  of  pulses 

•  optical  memory 

•  low-loss  bending  of  light 


The  I’nivcrsil)  »f  Nc\»  Mcxkh 


1d  strucutres:  Optical  fiber  gratings 


Fiber  Gratings 

Refractive  index  cf  core 
(germanium  doped 
silicate)  can  be  varied  by 
exposing  to  ultraviolet 
light  (laser  approx 2^ 

nm)  - 

- ◄ - 


Cladding 


Core 


Periodic  spacing  produces 
Bragg  grating  within  fibre.  The 
period  can  be  precisely 
determined  by  laser 


Ej-  (Z,  T)  (forwardmoving  envelope) 


Ej^  (Z,  T)  (backwardmoving  envelope) 


Equations  studied 


1-dim  Coupled  Mode  Equations 

=  -c^d^E^  +  iKE_  +  /r(|  E^  I"  +2 1  E_  f)E^ 
d,E_  =  c^d^E_  +iKE^+  iTCI  E_  f  +2  |  E^  f)E_ 


Gap  solitons  exist.  Velocity  proportional  to  amplitude  mismatch 
Between  forward  and  backward  envelopes. 


2D  structures:0  Waveguide  gratings 


Assumptions: 

•  dynamics  in  y  arrested  by  a  fixed  n{y)  profile 

•  xy-normal  incidence  of  pulses 

•  characteristic  length  scales  of  coupling, 
nonlinearity  and  diffraction  are  in  balance 


BARE  2D  WAVEGUIDE 


WAVEGUIDE  GRATING 


-  2D  NL  Schrodinger  equation 

-  collapse  phenomena:  point  blow-up 


-2D  CME 

-  no  collapse 

-  possibility  of  localization 


Ihc  I’nivcrsil)  of  Ntr.  Mexico 


Dispersion  relation  for  coupled  mode  equations 


Close,  but  outside  the  gap 
Well  approximated  by  the 
2D  NLSE  +  higher  order 
corrections. 

Collapse  arrest  shown  by 
Fibich,  Han,  A.A  (2002). 
But  dynamics  is  unstable 


Frequency  gap  region. 

We  will  study  dynamics  in  this 
regime. 


Governing  equations  :  2D  Coupled 

Mode  Equations 


lE^  -\-iKE_  +/r(|  Ej^  ^  +2  E_  P)-£' 


+ 


a,£  = 


c^d^E_  "^ifiE^  +/r(|  E_  +2  E^  \  )E_ 


,d,K,T  >(),  E+  \  [— L^,Z^]x[— Z^,Z^]x[0,oo)  — >  C. 


Summary  of  results  (details  to  follow 
on  next  slides) 


•  Found  stationary  solutions  of  the  governing 
equations  of  the  2d  structure  (see  next  slide) 

•  Obtained  conditions  for  bullet  propagation  in 
such  structures 

•  Obtained  conditions  for  light  trapping  at  a  defect 
by  a  resonance  mechanism  between  the  incident 
optical  bullets  and  defect  modes 

•  Derived  a  finite  dimensional  dynamical  system  to 
study  the  dynamics  inside  the  trap 


stationary  solutions  via  Newton’s  iteration 


If  E±{x^  2,  t)  =  S±{x^  z)  e  then 

ujE-\-  +  icgdzS-\-  +  d'^E-\-  +  t^E—  +  rd^^l"  +  2 \E— 1"")^+ 
ljE—  —  icgdzE—  +  9~E—  +  k,E-\-  +  rd^—l"  +  2 1£’-|-|  *')£’_ 


0, 

0. 


Solve  (1)  as  a  NL  eigenvalue  problem  for  (  a;, 


via  Newton’s  iteration. 


Need  one  more  eciuation: 


/oo  roo 

/  \S+\^+\e-\^dzdx  =  N. 

-OO  — oo 


Initial  guess:  E^\x, 

separable  wavefomi  E^\x,  z)  =  E±{z)G{x),  =  kcos{S) 

where  E±{z)  is  the  ID  gap  soliton  with  r?  =  0  (free  parameter  G  (0,  tt)) 

Substitute  and  integrate  in  z: 

G"  +  b{G^  -  G)  =  0,  h  =  2^{sm(S)-Sms(S)),  c5€{0,7r/2) 

() 


=>  G{x)  =  \/2  sex'h{VT)x) 


stationary  case  I 


Peak  amplitude  evolution 
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Nonexistence  of  minima  of  the  Hamiltonian 


Theorem  1.  The  Hcmultorium  funcJiojial  of  2D  CME  has  no  rninirna,  constramed  to 
total  po'U)ej\ 

Proof:  Suppose  z))  s.t.  H{E+,E_)  =  miusH^  where 

2 

•S’  =  {(Mi\z),f2(2\z))  S.t.  fi-2  :  ^  ^  II/A.II2  =  11^+112+  11^-11*1} 

k=l 

Consider 

Si  =  {{E^,E_)  :  E±  =  aE±{x/ p.  z/ ly)  and  a,//,//  >  (),a^////  =  1}. 
Clearly  S\  C  S  and  {E_^.,E_)  G  Si. 

Within  Si  the  Hamiltonian  H  =  H,.  =  —  .440^  —  .42  with 

.4i  —icg  I  EfdzE_  —  EfdzE^dxclz.  A2  =  n  I 

7r2  Jr 


.  A‘2  =  K  I  E_E\A  EfEj^dxelz, 

R2 


.43  —  f  \dj.E_^\‘^  +  \dj.E_\‘^dxdz,  .44  —  ^  f  \^+\^  +  \E_ 

Jm?  ^  v/r2 

.4i,  .42  G  M  and  .43,  .44  >  0.  The  only  C.P.  of  H,.  is 


(aC/y*)  = 


.4i 


.4f  ’2.4i.43 


a  fixed 


and  by  the  2nd  derivative  test  (a*,  lA)  is  a  saddle! 


□ 


Defects  in  1d  gratings 


■  Trapping  of  a  gap 
soliton  bullet  in  a 
defect,  (this  concept 
has  been  theoretically 
demonstrated  by 
Goodman,  Weinstein, 
Slusher  for  the  1-dim 
(fiber  Bragg  grating 
with  a  defect)  case). 


Position 


Ref:  R.  Goodman  et.al.,  JOSA  B  19,  1635  (July  2002) 


Governing  equations  : 
2D  CME 


+/3^2^+  -\-iK'(^x^z^E  -\-V (x,z)£^^  +/r(|  E_^  ^  +2  E 

d^E  =  3 ^E  +  /3  2^_  i z^E^  +  (^XyZ^E  +  /r(|  E  ^  +2 


advection  diffraction 


defect 


non-linearity 


AT,  r  >  0,  £,  :  [-4 , 4  ]  X  [-i^ ,  i  1 X  [0,  oc)  ^  c. 


2-D  version  of  resonant  trapping 


loft:  GS  modulus 

right:  defect  potoiitial  V[x,z) 


(-)  Bifurcation  curves  for  the  XL  defect  mod 
(*)  stationary  GS  with  u.’o  ~  0.96 


3  trapping  cases  (GS:  a)(v=0)  s  0.96) 

1 .  Trapping  into  two  defect  modes 

V  =  V,{x)T{z;9)  +  V^{z)T{x;l),  x-  =  -^l  +  /t'(tanh"(fe)-l), 

where  F[  =  2;ff^sech^(/?x),  ^2  -  ^  and  r(j^;c)  =  ^(taiih(>’  +  c)-tanh(>’-c)) 


2\  +  k  (tanh  (kz)  —  1) 


with  A:  =  0.1 8  and  =  0. 16 


2  linear  defect  modes :  co^  ~  0.963,  0.992 


3  trapping  cases  (GS:  (a)(v=0)  ~  0.96) 


2.  Trapping  into  one  defect  mode 


V  =  l  +  =  0.25, Z)  =  0.3 

^  1  linear  defect  mode  :  co^  -0.995 


E  (x  = 


3  trapping  cases  (GS:  (a)(v=0)  ~  0.96) 


3.  No  trapping 

V  =  V,{x)T{z\S)  +  V^{z)T{x\5\  k= -.j\  +  k\x?ir^\kz)-\\ 

where  =  2y^^sech^  {fix),  V^=l-  f  7  ^  1^2  ^  =  \  (tanh(j;  +  c)  -  tanh(  j;  -  c)) 

2\  +  k  (tanh  {kz)  - 1)  2 

with  k  =  0.85  and  13  =  \ 

^  5  linear  defect  modes :  -  -0.47,  0.26,  0.47,  0.68  and  0.87  ^  ALL  FAR  FROM  RESONANCE 


E^{x  =  X),z,t)\  l«^(OI 


Finite  dimensional  approximation  of  the  trapped  dynamics 

For  trapped  solutions  with  small  amplitude: 


E+{x,z,  f) 
eJ{x,z,  t) 


N 


where  ’>  A’  =  1, .  • .  are  the  defect  modes. 

Substituting  into  2D  CME  with  the  defect  potentials 


where  [NL)±  =  {\E± 
Due  to  orthogonality 


N 


k=l 

+  2EX-)E±. 


ia'X)  +  r  /  (NL)+i>*^^  +  {NL)-xl)*_^dxdz  =  0  , 


A-  =  1. . 


*  *  •> 


N. 


For 


^  Finite  dimensional  approximation  of  the  trapped  dynamics 


after  transformation  ni,2(f)  =  where  Auo’ 

Ao/'  .  9  _  1 1  ~  ~  ~  9  ~  * 

“Oi  +  ;i|n2|*”ni  + 


la 


1 - ^<^^1  +  <^i 


iCL2  ^2  “1“  ^2  *^^2  “1“  ,'^1^1  ^^2  “1“  ^^1^2 


0 

0. 


Conclusions  (Part  I) 


■  2d  nonlinear  photonic  structures  show 
promising  properties  for  quasi-stable 
propagation  of  “slow”  light-bullets 

■  With  the  addition  of  defects,  we  presented 
examples  of  resonant  trapping 

■  Research  in  line  with  other  interesting  schemes 
to  slow  down  light  for  eventually  having  all 
optical  logic  devices  (eg.  buffers) 


Part  II:  Dynamics  in  Bose- 
Einstein  Condensates 


statement  of  the  problem: 

In  the  study  of  the  dynamics  of  matter  waves  for 
Bose-Einstein  condensates,  it  is  of  great  importance  to 
understand  the  scenarios  under  which  solutions  such  as 
necklaces,  multi-poles,  and  vortices  will  exist  and  persist. 


Part  2:  Bose  Einstein  Condensates 
(summary  of  most  important  results) 


■  Shown  that,  under  suitable  assumptions,  that  N-solitons 
are  stable  for  a  large  class  of  integrable  partial 
differential  equations  including  the  equation  that  govern 


BEC 


■  Shown  how  the  presence  of  an  optical  lattice  along  with 
the  magnetic  trap  can  influence  the  dynamics  of  matter 
waves  in  Bose-Einstein  condensates 


■  Illustrated  the  manner  in  which  increasing  the  number  of 
potential  wells  in  photorefractive  media  effects  the 
dynamics  associated  with  solitary  waves 

■  Begun  to  develop  some  "rules-of-thumb"  regarding  the 
existence  and  stability  of  waves  in  two-dimensional 
Bose-Einstein  condensates 
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